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We study the zero-temperature spin fluctuations of a two-dimensional itinerant-electron system 
with an incommensurate magnetic ground state described by a single-band Hubbard Hamiltonian. 
We introduce the (broken-symmetry) magnetic phase at the mean-field (Hartree-Fock) level through 
a spiral spin configuration with characteristic wave vector Q different in general from the antiferro- 
magnetic wave vector Qaf, and consider spin fiuctuations over and above it within the electronic 
random-phase (RPA) approximation. We obtain a closed system of equations for the generalized 
wave vector and frequency dependent susceptibilities, which are equivalent to the ones reported 
recently by Brenig. We obtain, in addition, analytic results for the spin-wave dispersion relation in 
the strong-coupling limit of the Hubbard Hamiltonian and find that at finite doping the spin-wave 
dispersion relation has a hybrid form between that associated with the (localized) Heisenberg model 
and that associated with the (long-range) RKKY exchange interaction. We also find an instability 
of the spin-wave spectrum in a finite region about the center of the BriUouin zone, which signals a 
physical instability toward a different spin- or, possibly, charge-ordered phase, as, for example, the 
stripe structures observed in the high-Tc materials. We expect, however, on physical grounds that 
for wave vectors external to this region the spin-wave spectrum that we have determined should 
survive consideration of more sophisticated mean-field solutions. 

I. INTRODUCTION 

Spin fluctuations about incommensurate static magnetic configurations represent an interesting problem, especiallst 
in the view of the accumulating experimental evidence on the parent compounds of the high-Tc cuprate materials.EJ 
The dynamic structure factor S'(q, oj) at wave vector q and frequency w, as measured by inelastic neutron scattering, 
shows noticeable peaks over the background when q = Q = Qaf + A Q (|A Q| ^ |QAF|)nfor increasing dopant 
concentration, at finite albeit small values of ui (with a; <C J, J being the exchange coupling) .□ Further experimental 
studies for larger values of |q — Qaf| have also detected the presence of well-defined spin- wave- type excitations close 
to the boundary of the antiferromagnetic (AF) BriUouin zone (BZ).B The fact that spin- wave excitations with small 
wavelength (i.e., comparable to lattice spacing) can be detected even in the absence of long-range magnetic ordei', has 
actually been well established experimentally since the eighties also in more conventional magnetic materials.cl This 
state of affairs has, iaijurn, prompted a number of theoretical studies on dynamical excitations about incommensurate 
spin configurations .ETEI 

From the theoretical point of view, inclusion of incommensurate spin configurations in cluster mean-field calculations 
would require one to consider cluster sizes at least as big as the spatial extension of the incommensurability. In this 
way, complex incommensurate patterns could as well be included in the calculation, at the price of considerable 
numerical effort and without full analytic control on the results. Alternatively, one may set up calculations for an 
infinite system with necessarily simpler incommensurate spin configurations, with the advantage, however, of obtaining 
analytic results (at least) in some limits, which in turn may admit a simple physical interpretation. In this respect, the 
' spiral spin configuration appears to be the only one for which analytic calculations can be performed, in the sense that 
the constituent equations can be brought to a closed form which is manageable for controlled analytic approximations. 
Indeed previous analytic calculations have considered stoti«|4ong-range background spiral spin configurations, on top 
of which dynamical spin excitations have been considered .BlZI 

The intrinsic simplicity of the background spiral spin configuration and the neglect of its dwiamics have resulted, 
however, into an instability of the spin-wave spectrum in a limited region of the BriUouin zoncQ This outcome could, 
in principle, make the resulting spin-wave spectrum altogether unreliable, as the underlying spin configuration (over 
and above which spin-wave excitations are constructed) could be much more complex than a spiral one and even 
possess a dynamics of its own. In addition, a static or dynamic charge modulation can be present as well, like 
the stripe structures observed in some underdoped cupratesJj One expects, however, on physical grounds the short- 
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wavelength spin-wave excitations, obtained on top of a spiral configuration, to preserve their dispersion relation even 
when considering more complex long-range underlying structures. 

Within this framework, we have pursued the analytic calculation with an underlying spiral spin configuration for an 
itinerant electronic system described by a single-band Hubbard Hamiltonian. The spin (and charge) dynamics have 
bee described within the electronic random-phase (RPA) approximation for the dynamical susceptibilities, based on 
a broken-symmetry Hartree-Fock (HF) mean-field solution, along the lines of RefE. In particular, wc have considered 
the large U/t limit in detail (where t is the nearest-neighbor hopping matrix element and U is the local on-site 
repulsion in the Hubbard Hamiltonian) , for which the spin- wave dispersion relation (as obtained from the poles of the 
dynamical susceptibilities) can be expressed in analytic form by systematically expanding in the small parameter t/U . 
In this way, we have obtained an analytic expression for the spin- wave dispersion relation valid in the limit t/U <^ 1 
[cf. Eq. (|5^) below], showing a novel characteristic structure, namely, a hyb«d form between the dispersion relations 
obtained within the nearest-neighbor jHeisenberg model for localized spinal^ and the long-range RKKY interaction 
mediated by the conduction electrons O 

This analytic form of the dispersion relation could admittedly not have been guessed a priori, by fitting the 
dispersion relation obtained numerically (from the location of the poles of the dynamical susceptibilities) with a 
Heisenberg model extending in principle to a large albeit finite number of neighbors. Neither, this analytic form cau 
be simply reduced to a nearest-neighbor Heisenberg dispersion relation with a doping-dependent exchange integral.E3 
Rather, the characteristic long-range RKKY contribution would require fitting to a Heisenberg model with an infini|te 
number of neighbors. This would contradict the spirit with which the Heisenberg model was introduced to start with,Ej 
namely, as a fitting model that makes physical sense when the interactions extend to a limited number of neighbors 
only. Note that this situation contrasts that found at half- filling of the Hubbard band (i.e., in the absence of doping), 
where the antiferromagnetic (AF) spin-wave spectrum can be nicely fitted by a Heisenberg model extending at most 
to a few neighborsJlj Our results also show that the magnitude of the overall exchange integral, which characterizes 
the spin-wave spectrum, decreases with increasing doping and vanishes when the transition to a ferromagnetic case 
occurs at the mean-field level. 

The plan of the paper is as follows. Section 2 obtains the implicit form of the spin- wave dispersion relation within 
the HF-RPA approximations, by solving for the dynamical susceptibilities of the itinerant electron system in the 
presence of an incommensurate spin spiral ground state. Section 3 focuses on the small t/U expansion of the results 
of Section 2, which requires a careful analysis of the doping dependence of the relevant quantities. Section 4 discusses 
the main results of this paper and Section 5 gives our conclusions. For the sake of completeness, we report in the 
Appendices details of the analytic calculations, and adapt know results for the Heisenberg and RKKY spin-wave 
spectra to the present context. 



II. DYNAMICAL SUSCEPTIBILITIES WITHIN THE ITINERANT-ELECTRON RPA 
APPROXIMATION WITH AN INCOMMENSURATE SPIN-SPIRAL GROUND STATE 

In this Section, we give the derivation of the spin-wave dispersion relation for a two-dimensional itinerant-electron 
system in the presence of an incommensurate spiral magnetic structure with a generic characteristic wave vector 
Q, within the electronic RPA approximation. Although our results for the dispersion relation coincide with those 
previously given by Brenig in Ref-Q, we provide here in addition the expression of the correlation functions which 
can be relevant for a direct comparison with neutron scattering experiments. Further details of the calculation are 
reported in Appendix A. 

We emphasize that the finding of a closed-form expression for the correlation functions for any characteristic wave 
vector Q (incommensurate with the lattice spacing) and not just for the (commensurate) antiferromagnetic wave 
vector Qyi_F is altogether a nontrivial result, being intrinsically related to the peculiar pattern of the spiral magnetic 
solution for the ground state. 

We begin by considering the generalized correlation function at zero temperature in the broken-symmetry phase: 

X^^,{rt,r't') = -t{T[S^{r,t)S,{r',t')])+t{S^{r,t)) (5.(r',i')) (1) 

where the average (• • •) is taken over the ground state, T stands here for Wick's time-ordering operator, and S'^(r) is 
given by (we set h = 1 throughout) 
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(2) 



a, 13 



In these expressions, ^x^v = {0,x,y,z), is a Pauli matrix (with cr° equal to the 2x2 identity matrix), and a,/? 
are spin labels. Note that 25o(r) coincides with the density operator, which couples with the spin operator in the 
presence of an incommensurate spiral magnetic structure. 

For the simple band we are considering, the field operator in Eq.(|^) can be represented in the form 



(3) 



where 0(r) is the atomic (Wannier) orbital associated with the simple band, R,; is the lattice vector locating site i, 
and is a destruction operator. Time evolution in Eq. (m) is governed by the Heisenberg representation: 



■ipair, t) 



iHt 



V'a(r) e" 



iHt 



where for H we take the simple-band two-dimensionaUHubbard Hamiltonian. 

In terms of the two-particle correlation function L,Ej the generalized correlation function (Jl^) takes the form 



A'^,.(rt,r'i') = E E <P L{rtf3yt'/3';rt+ayt'+a') 

a, 13 a' ,13' 



with i+ = 1 + 7] (r/ = 0+), where L satisfies the Bethe-Salpeter equation: 

L(1,2;1',2') = G(1,2')G(2,1') 

+ J didAdb d6G(l,3)G(4, l')S(3,5;4,6)L(6,2;5,2') 

(1, 2, • • • signifying the set of space, spin, and time variables). In the above expression, 

G(l,2) = -z(r [^(l)^n2)]) 



(4) 



(5) 



(6) 



(7) 



is the single-particle Green's function and the kernel ^ represents an effective two-particle interaction. In particular, 
within the RPA approximation we are adopting, the kernel S takes the form: 



S(3, 5; 4, 6) = -ivq i5(3, 4) (5(5, 6+) (5(x3, xe) (5(a3, -ae) 
-t- Wo (5(3, 6) (5(4, 5) 5{x^ , X4)5{a^, —ai) 



(8) 



with the notation x = (r, t) and where the constant vq can be related to the parameter U of the Hubbard Hamiltonian 
as follows: 

U^voJdv\cl>{r)\^. (9) 

Entering Eqs.(^ and (^) into Eq.(^), we obtain for the generalized correlation function within the RPA approxi- 
mation: 



'^/i,l/('^5 ) 



'iE E <,p<',p'G{xP,x'a')Gix'(3',xa) 

a, 13 a' ,13' 
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H)'EE <P<'.P' [d3Gix(3,xsa:,)Gi 



X3a3,xa) 



a, 13 a', 13 



T(-*)'EE <P<':P' fd3G{xP,X3a3)G{. 

n R n' Rl 



a;3Q3,a;Q!) 



a, 13 a', 13' 



X L(x3a3, a:'/?'; ^3 a3, x'^a') 



(10) 



with a = —a. This is apparently not a closed-form equation for X itself. By the manipulations reported in Appendix 
A, however, Eq. ( pO| ) can be cast in the form of a coupled set of equations for the matrix components of the correlation 
function, as follows: 
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X^^,{x, x') = A-W (a;, a;') + 2 «o 51 j ^x" -^^"i {x, x") t^,. y X,,,,{x'\ x') 



(11) 



where we have introduced the tensor 



/I 

0-100 
0-10 

Vo 0-1 



(12) 



To solve Eq. (jTTI), the explicit form of the non-interacting counterpart A"*^*^-* of X is required. To this end, the 
ground-state average in Eq. (|^) is evaluated as shown in Appendix A, in terms of the eigenvalues (e^) and eigenvectors 
(W^,r) of the mean-field Hubbard Hamiltonian (see also Section 3), yielding for the time Fourier transform of X'^^'^ 
the expression 



BZ 

<Y(«)(r,r';c.) = e^^'^-'^'H^-^.) |0(r-R,)P |0(r' - R,)|^ 

i,j k,k' 

EE ^M.M'(f^or.,.'(f^,)F/t;(k',k)F;::,(k,k') 

X !F^y (k, k', Lo) 



(13) 



where BZ stands for the two-dimensional Brillouin zone, M is the number of lattice sites, and the quantities T, _F, 
and T are defined in Appendix A (17,; standing for the angles defining the local spin quantization axis). We consider 
further the space Fourier transform 



A',b(q,q';co) - j^J dvdv' e-'^ '' Xab{v,v'-u:)e'^-'' 



5(q) S*{c{)Xab{ci,c{-Lo) 



where Vo is the volume of the elementary crystal cell and 



5(q) 



(ire 



(14) 



(15) 



is a form facto r (which can be set equal to unity for all practical purposes). Applying a suitable unitary transformation 
[cf. Eq. ( A23| )] which renders the matrix T{VLi) of Eq. ( p^ ) diagonal, one gets for its lattice Fourier transform: 



Tafc(k) = E ^"^'^^ ^-"{n^ 



/A(k) 

A(k-Q) 

A(k) 

V A(k + Q) 



(16) 



A(k) being the lattice Kronecker delta function. In the new basis, we thus obtain for the non-interacting correlation 
function the expression: 

BZ 

(q. q'; '^) = 777 E E E q) ^'^'^(k' - k + q') 



kk' a' ,b' r,r' 
7b' 



X F;, ,(k', k) F,':^, (k, k') Try (k, k', c^) 



(17) 



where the overbar denotes matrices transformed according to the above unitary transformation [cf. Eqs. ( A19 ) and 
([A25|)]. In this way, the integral equation (O) reduces to the form: 



A',b(q,q»=<Y(°)(q,q' 



2t^E E ^aa' (q' q"; ^) -^ff (q"' q'; 

a'h' q" 



(18) 
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with e given by Eq. (A22). This equation can be solved by the methods of Appendix A, yielding the closed- form 
expression: 



Xab{ci + 7aQ, q'; u;) = J2[l + X(q, lu)]-^, 



Xl°l (q + 7a'Q; co\Q) A(q - q' - 7faQ) 



(19) 



where 1 is the 4x4 unit matrix, the matrix X(q, w) is defined by Eq. ( A30| ), and with the notation 7a = for 
a = 0,2, 7a = —1 for a = 1, and 7a = 1 for a = 3. Although still expressed in the transformed basis, Eq. (^9|) is 
the desired expression for the Fourier transform of the generalized correlation function, which holds within the RPA 
approximation for any value of Q. 

The spin-wave dispersion relation can eventually be obtained by searching for the zeros of the inverse matrix on 
the right-hand side of Eq. (fol), which is equivalent to imposing the condition: 



det[l-|-X(q,w)] = 0. (20) 

It can be verified that the condition ( |20| ) can be mapped onto the result reported in RcfE, where the dispersion relation 
has then been obtained numerically for chosen values of Q. In the present paper we proceed instead to deriving the 
analytic form of the dispersion relation for small values of the parameter t/U of the Hubbard Hamiltonian. 

To this end, it is convenient to rewrite first the matrix X in Eq. (|2^) in a more conventional basis identified by the 
labels (0, +, — , z), with cr* = (cr^ ± z(t*')/\/2- The matrix 1 + X(q, uj) is then transformed into: 



/ 



M{q,uj) = 1 + 2C/ 



■vO.Q/ \ 



V 



-X-''i<i,cj) X- 



X°^-{(l,uj) 



-A'o^^"(q,c.) 



Xn 



'(q^w) 



A'°'+(q,c.) 



Xq 

Xn 



'+(q,^) 
'+(q,t^) 



Xq '^(q,cj) 



(21) 



where now 



1 

AN 



BZ 



^ ^ F,",^,(k - q, k) F^^, (k, k - q) T^y (k, k - q, c.) 
k 



(22) 



and 



E 



(k)af_^,t^e 



'(k') 



(23) 



with (7 



+ 



^/2 ^ Q Q J a = \/2 ^ ^ g J . Note that two columns in the expression (|2l|) appear interchanged 
with respect to the original order, owing to the presence in the final basis of the tensor 




(24) 



in the place of e given by Eq. ( A22| ). 

We pass now to perform the small t/U expansion of the spin- wave dispersion relation obtained from the condition 
detM(q,u;) = 0. 
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III. EXPANSION IN THE SMALL PARAMETER t/U 



In this Section, we obtain explicitly the spin wave dispersion relation to second order in the small parameter 
t/U, from the general condition detM(q, lj) — obtained in Section 2 within the RPA approximation for the zero- 
temperature broken-symmetry phase with a generic incommensurate wave vector Q. To this end, we will preliminary 
expand the self-consistency parameters of the mean-field Hamiltonian, as well as the equations they satisfy, at the 
relevant order in t/U] we will then expand the matrix elements of the matrix Af(q, cj) defined by Eq. (|2l| ) at the 
relevant order in t/U . 



A. Mean-field equations 



The mean-field equations for a single-bapd Hubbard Hamiltonian in the presence of an incommensurate spiral 
spin structure have been discussed in Ref.tS. Introducing a local set of spin quantization axis, with the z axis 
transformed locally into the axis specified by the spherical angles 51^ = (0^ = Q • R^, (^^ = 0) at site i, one transforms 
the destruction operators Cia according to Eq. ( |A6| ) and performs the (Hartree-Fock) mean-field decoupling of the 
Hubbard Hamiltonian, yielding 



BZ 



with ^, ^' — (+, — ) and where H is the 2x2 matrix 
H(k) = 



eo - fi + tTeik) + U {mi - m2) -itTo(k) 

itTo(k) eo ~ fi + tTeik) + U {mi + m2) 



(25) 



(26) 



In the above expressions, eo is the site energy, /i the chemical potential, mi and m2|-i:epresent the occupation number 
and magnetization along the local quantization axis, respectively, and Tg/o(k) readll3 



TJk) 



cos fc-r cos 



cos fc„ cos 



Qy 



(27) 



TJk) 



sinfca; sin 



sin ky sin 



2 J " V 2 

The eigenvalues and eigenvectors of the matrix (M) are thus given by: 



(r = 1,2) and 



e,(k) = eo-^l + tT,{k) + Umi + {-IfJU^ml + f2To{kY 



(28) 



(29) 



M^i(k) = 



iVi(k) 



(30) 



W2(k) 



7V2(k) 



'i+^/i+(ir^^^Mk)^ 



(31) 



where N^ik) stands for the normalization factor. 

The parameters mi, m2, and Q of the mean- field Hamiltonian ( p5| ) are obtained, as usual, by minimizing the 
average value of the Hamiltonian with respect to the parameters themselves. One obtains: 
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1 

"^1 = /f(er(k)) 



l + (5 



(32) 



k r 



where /i?(e) is the (zero-temperature) Fermi function and 5 is the doping parameter^ 



BZ 



k C r 



and 



BZ 



^Q^«.«'(k) ^<5(k)W-e,r(k)/F(e.(k)) = 

k 5,4' r 



(33) 



(34) 



In the following, w&shall restrict to the diagonal solution Q = Q(l, 1), since it is known to be favored for sufficiently 
small values of t/U^3 Accordingly, at the order we are considering of the small parameter t/U we expand formally: 



4")(k)-M(") 

2 



e^(k) = U 

+ (^)'(4^Hk)-^(^)) + .. 



(eW(k)-MW) 



as well as 



(0) , 

7712 — 'm-2 + 



,(2) 



(35) 



(36) 



while, at the relevant order we can take 



„■ t T„(k) 



(37) 



1 /_i±IkM\ 

W2{k) = -^i *^2™<«) (38) 



where 



Nriky 




(39) 



is independent from i\ Note that it is sufficient to retain the lowest-order term 171^"^ in the above equations. The 
parameter mi, on the other hand, is given by Eq. (j3^) and is thus formally independent from t/U (we anticipate, 
however, that the doping parameter 6 will turn out to be at most of the order t/U for our expansions to be internally 
consistent). 

The coefficients e["^(k) of Eq. (H) with 77 = 0, 1, 2, • • • can be readily obtained from the expression (^ ) (where we 
may set eg = for simplicity) in terms of the 7712"" (for given Q). The value of /i'-"-' can also be readily obtained in 
terms of e[*'^ (cf. Appendix B). The remaining coefficients of the expansions ( |35| ) and ( |36| ) can further be determined 
by solving the coupled equations for the self-consistency parameters according to the methods developed in the 
Appendices B and C. In particular, for S > we obtain: 

777(°).i(l-^) 



777^^^ = 



-4sin2(Q/2) + 0((5) , (40) 
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(0) 




(1) 

— 1712 


^(^) 


(2) 



= cos(g/2)[4-47r(5 + C'((52)] 

= + O(<5), (41) 



and 



e(°)(k) = i[(l + 5) + (-l)'-(l-<5)] 



2 

e«(k)=re(k) 

e?Hk) = i-iy [T2(k) - 4 sin2(Q/2)] + 0(<5) , (42) 
where for the spiral configuration we are considering Q is determined by the condition 

cos(Q/2) = ^+0 {t/U) (43) 

with 5 < 2{t/U) + 0{{t/U)'^), as anticipated. (The other allowed solution sm{Q/2) — describes instead the 
ferromagnetic case.) As discussed in Appendix B, the above results have been obtained with the further assumption 
that 6 is small enough but not infinitesimal, i.e., 5 satisfies the condition {t/UY « The case (5 = 0, on the other 
hand, can be considered separately. Note that Eq. ( |4^ ) implies that in the spiral phase 6 is at most of the order t/U . 
This property has to be taken into account to get a consistent expansion up to the desired order in t/U . 



B. Susceptibilities and spin-wave dispersion 



Before performing the t/U expansion of the matrix elements of the non-interacting susceptibility tensor ( p^ ) to get 
the spin-wave dispersion relation, it is convenient to exploit some symmetry properties that reduce the number of 
matrix elements to be considered. Specifically, from the property 



J^r,r'(k,k- q,w) = Tr',rO^~ q, k, -w) 



it follows that: 



By direct inspection it can also be verified that: 



"(q,^) = -Xo'''ici,uj) , 
"(q,^) • 

In this way, the matrix (|21j) acquires the simplified form: 

i5(q, iv 

M(q,w) 



/ l-a(q,tj) 

ib{q^,uj) l + d{ci,ijj) 
-i6(-q,-cj) e(q,u;) 



e(q,w) 
1 + d(-q, -uj) 

«/(-q, 



where we have set 



a{q,uj) 


-2[/A'°'°(q,o.) 




6(q,tj) 


= -2zC/A'o"-(q, 




c(q,w) 


^2UX;^'\q,iu) 




d(q,tj) 


= 2UX+'-ici,Lu) 




e(q,tj) 


= 2UX+'+i<i,u;) 




/(q,t^) 


^-2iUX+^\ci, 


c) 




^2UX^^'{<l,u;) 





(44) 
(45) 



c(q,tj) 
i/(q,t^) 

-«/(-q, 
1 + g(q,t^) 



(46) 



(47) 



(48) 



8 



Entering then the expansions (|3^)-(|3^) for the eigenvectors of the mean- field Hamiltonian (with mlp given by 
Eq. (^)) into the definition (^3|), we obtain the expressions for the relevant matrix elements of the non- interacting 
susceptibility tensor (^) reported in Appendix B at the order in t/U we are considering. Utilizing further the method 
developed in Appendix C to perform the k summation when the doping parameter 5 is small, we obtain eventually 
the following expressions for the matrix elements 



as well as 



a(q, w 


a(q) 


6(q,w 


= fe(q) 


c(q,w 


-c(q) 


e(q,tJ 


= e(q) 




= /(q) 


5(q,w 


= 5(q) 




1 - a) + 




1 - a) + 


+ 2 




- 2 





1 



5U 



2 cos((3/2) [cosqx + cos qy — 2) t 
- sin((5/2) {sinqx + sing^) 



O (t/[7) , 



\/2 cos(Q/2) (cosg^ + cosg.y — 2) 



(5 + O (vc/) 



-a(q,cj 



O (t/C/) , 



: 6(q,w) 
a(q,tj) 



sin (Q/2) (cos(7j: + cos qy 



o[it/uy 



(49) 



— j 4 cos^(Q/2) (cosg^ 



COS( 



2) + 8 - sin2(g/2) 



+ COS qy — 2) 



sin^((5/2) (sinQj; 



iqy) 



Uj cos{Q / 2) {cosqx + cos qy — 2) 



O 



(mf) 



(50) 



where uj = ui/U will turn out to be of order t'^ /U'^ at the spin- wave poles. 

To obtain the above expressions, we have considered only the real part of the functions ( A16| ). This is definitely 
possible for every pairs of bands when we restrict to values of q such that |q| 3> kp, where kp (by our definition) is 
the maximum value of the function k{4>) introduced in Appendix C, which coincides with the Fermi momentum for 
small S. Since we have shown in the same Appendix that k{(jji) ~ ^/S ~ ^Jt/U, taking the q ^ limit implies letting 
t/U to vanish before q. 

Note also that, although the expressions ( |4£) a nd ( ^0| ) have been calculated at different orders in t/U, the ensuing 
expression for the determinant of the matrix (4^) is obtained consistently at the fourth order in t/U, as required for 
the frequency of the spin- wave mode to be of the order ot t^/U. In fact, by writing the determinant explicitly we 
obtain: 

[(1 - a(q)) (1 + a(q)) + a^iq) + O {t/U)] 

X \-u + a(q)) + a(q)) - e2(q) + O ({t/Uf) 

+2 [1 - a(q) + O {t/U)] [fe2(q) + O ((t/C/)')] [e(q) - «(q) + O ({t/Uf 
-2 [1 + a(q) + O {t/U)] [52(q) + O ({t/U) 



^A[a{q) + 0{t/U)] 



O 



[{t/U)') 



e(,q 



e(q) - a{q) + O [{t/U) 
)~a{q)+0({t/Uf) 



(51) 



Note that the last three terms on the left-hand side add up to zero at the fourth order in t/U we are considering. We 
are thus left with the expression 



(^2_a2(q) + e2(q) + 0((t/[/)^) 



= 0, 



(52) 



yielding 
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(53) 



Taking into account the expressions for a{q) and e(q) given above, we get eventually the desired spin-wave dispersion 
relation, in the form: 

w^(q) — 16 — cos Q (cos Qx + cos qy — 2) [cos qx + cos qy ~ 2 cos Q] 



cos((5/2) (cosQj; + COS qy 



sin^(Q/2) (sin^a; + sin^j^)^ 



cos(Q/2) (cosgj; + cos qy — 2) 



X [cos^(Q/2) {cosqx + cos qy - 2) + 2 sm^{Q/2)] 



cos((5/2) (cos (7a; + COS qy — 2) — 



siv?{Q/2) (siuQa; + sin9y)2 
cos((5/2) (cos + cos qy — 2) 



(54) 



To obtain the physical dispersion relation, however, there remains to enter in the above expression the relation among 
Q, 6, and t/U as given by the mean-field condition (E^), as discussed in the next Section. 



IV. RESULTS AND DISCUSSION 



In this Section, we discuss the physical consequences of the dispersion relation (|5J) for spin-wave excitations over 
an incommensurate (diagonal) spiral magnetic configuration of a two-dimensional Hubbard Hamiltonian. 

We have already remarked that the dispersion relation ( |5^ ) is not yet in its final form, since the connection among 
Q, 6, and t/U needs still to be specified. Before considering the general case, however, it is interesting to recover 
from Eq. (m) the spin-wave dispersion relations corresponding to the limiting cases of an antiferromagnet and of a 
ferromagnet. 



When 6 = 0, Eq. ( |B14D yields cos((5/2) = 0, that is Q = Qaf = (7'',7r). In this case we obtain from Eq. ( |54| ) (by 
setting 6 = identically therein): 



w(q) 



jfff'^^ A- {cos qx + cos qyf 



(55) 



with Jq^^^"* = 4i? /U . This result coincides with the spin-wave dispersion relation of a two-dimensional Heisenberg 
antiferromagnet at leading order in t/f/.tZI 
When Q = {2tt, 27r) and 6 arbitrary, we obtain instead from Eq. ( j^) 



uj{q) 



1 - 



5U_ 
'2t 



(2 - cosqx ~ cos qy) , 



(56) 



which coincides with the spin- wave dispersion relation of a nearest-neighbor Heisenberg ferromagnet with J^g- = 

At'^/U [1 ~ 6U/{2t)]. Note that J^'g'' is negative since the (mean-field) ferromagnetic solution is actually stable when 
6 > 2t/U 



In the general case of a (diagonal) spiral spin configuration, the relation among Q, 6, and t/U is given by Eq. (B14), 
that is, 6 = —{2t/U) cos Q/2 + 0{{t/U)^) at the leading order in t/U. Eliminating Q in favor of 6 and t/U via this 
relation in Eq. (p4), we obtain eventually the following expression: 



t^(q) = Jeff 



2- 



(smqx 



2 — cos qx - 

where now the effective exchange integral is given by 



COS( 



— (cosq^r + cos qyY 



(57) 
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with Sc = 2t/U. Equation ( p7|) constitutes the main result of this paper. Note that when S reaches the critical value 
<5cj JcS vanishes. Past this value, the spiral solution evolves into the ferromagnetic solution, which becomes the stable 
solution (at the mean- field level). 

The real and imaginary parts of i-u{c{) (in units of Jeff), as obtained from the analytic expression (|57|), are plotted 
in Figs. 1(a) and 1(b), respectively, over the two-dimensional Brillouin zone. 




diagonal-spiral spin configuration background. Note that the Q-dependence is contained in J^b only. 



It is clear from the expression ( p7[ ) that uj{q) is either real or purely imaginary, so that, when its real part 
is nonvanishing, its imaginary part vanishes identically, and viceversa. Actually, this is true in the region of the 
Brillouin zone where our expansion holds, namely, for |q| ^ kp [cf. the discussion below Eq. (pjQ^]. For |q| < kp, 
the spin-wave dispersion acquires a damping due to the mixing with the particle-hole continuuniEl An exception is 
represented by the line Qx = Qy, along which the real and imaginary parts vanish simultaneously. The softening of the 
dispersion relation along the whole line = Qy and not only when q = and q = Q (as one would expect on general 
grounds in the presence of a spiral spin configurationQ) corresponds to the fact that w(q) given by Eq. (|7|) can be 
cast in the form a;(q) = /(Q) ffQ(q); as a;(q = Q) = invariably, it follows that 5Q(q = Q) = since gqiq) does not 
depend on |Q|. This implies that ggiqx = Qy) =0 identically. 

One expects this result to be modified, however, at higher order in t/U. In this respect, it is interesting to compare 
with the results obtained by BrenigjZl by solving numerically the condition ( po|) directly, without performing the 
expansion in the small parameter t/U. One sees, in particular, from Fig. 6 of Brenig's paper (obtained for 6 = 0.075) 
that ti-'(q) remains indeed finite along = qy already when t/U = 0.1. 

Figure 2 shows the region of the Brillouin zone (shaded area) where w(q) is overdamped (purely imaginary) (with 
the region about q = excluded according to the argument given|in Section 3). A finite region of the Brillouin zone 
where the dispersion relation is overdamped is also reported in Ref.Ll, even though a direct comparison with our results 
is not possible owing to the different ranges of the parameter t/U explored. This overdamping signals an instability of 
the system (due to the merging into the particle-hole continuum) toward a different ground state, reflecting possibly 
a more complicated long-range spin (and charge) structure than the spin-spiral one considered in the present paper. 
Nonetheless, we expect on physical grounds that close to the boundary of the Brillouin zone (where overdamping of 
spin waves does not occur in our calculation) the spin-wave spectrum with small wavelength obtained by our approach 
would survive the inclusion of more complicated long-range spin structures. 

It is also interesting to compare the form of the dispersion relation (|5^) with the spin-wave dispersion relation 
obtained with the Heisenberg model including second and third nearest-neighbor couplings, for the same value of the 
incommensurability wave vector Q (cf. Appendix D). This comparison is shown in Fig. 3 along the symmetry lines of 
the Brillouin zone and evidences marked differences between the two dispersion relations. 

Returning to the dispersion relation (5^) for (5 ^ 0, we emphasize that its functional form could not be obtained 
from the dispersion relation (pq) valid when (5 = 0, by simply modifying the numerical values of the exchange integral 
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n 



-n 








n 



-n " q 

X 

FIG. 2. Region of the two-dimensional Brillouin zone where the spin- wave dispersion is overdamped (shaded area). Note 
that this region does not depend on Q, and that the region about the center of the Brillouin zone has been excluded by 
construction. 



therein. Nor, it would be sufficient to include a finite number of exchange integrals in the Heisenberg model to account 
for the finite doping, owing to the presence of an RKKY-type term in Eq. (|5^) which contains trigonometric functions 
in the denominator rather than in the numerator only (cf. Appendix E). Recall, in fact, that for practical purposes 
the Heisenberg model can be regarded as a fitting model which could, in principle, reproduce any type of spin-jsiave 
dispersion relation, provided a sufficiently large number of terms associated with distant neighbors were included.Ej In 
addition, it appears fair to say that it would have been certainly difficult to guess a priori the functional form of the 
dispersion relation (|5^), by fitting the numerical results for the itinerant model onto the Heisenberg model extended 
to a large number of neighbors. 

It is also interesting to comment on the effective exchange integral (|5|) retaining the t'^ /U dependence of the nearest- 
neighbor Heisenberg model, even when couplings between far apart neighbors are considered. From a perturbative 
point of view, when U » t the magnetic interaction between a given lattice site and far apart neighbors is provided 
by the mobility of the holes in the doped configurations. In this respect, one should consider all possible configurations 
with empty sites distributed at random over the lattice sites, the mobility of the holes then resulting by diagonalizing 
the Hamiltonian in this basis. As a consequence, the magnetic exchange coupling turns out to be proportional to 
t^/J7 at leading order, even for coupling between sites at arbitrary distances. 

As anticipated in the Introduction, the form (^7|) is somewhat hybrid between the one obtained with a nearest- 
neighbor Heisenberg model (cf. Appendix D) and the long-range RKKY interaction mediated by the conduction 
electrons (cf. Appendix E). These two contributions to Eq. (|5^) cannot be separated in a clear and unambiguous way. 
However, it is possible to trace their origin by considering the expressions ( B15| ) and (B16) for the relevant matrix 
elements of the susceptibility. 

If one could set "by hand" J-22 = in those expressions, thus keeping only the interband terms, one would in fact 
obtain for the dispersion relation: 



•'(q) = ^off [4 - (cosq^ + cos qy f] ' 



(59) 



which corresponds to a nearest-neighbor Heisenberg antiferromagnet with exchange coupling given by Eq. (|5q). Sim- 
ilarly, if one could set J^2i = J^i2 ~ "by hand" , thus keeping only the intraband terms, one would instead obtain an 
expression of the RKKY-type: 



(sing^ -I- sing;^)^ 
2 — cos qx — cos qy 



(60) 



with Jeff still given by Eq. (p8|). Note that in both cases the spectrum would be real and no overdamping would occur. 
Our general result (p^) can then be cast in the following appealing form: 



t^(q) - {[^(q-O) - c^(q)my]' - J(q)'} 



1/2 



(61) 
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2.0 

Re a)(q) 
J 

1.0 




(0,0) (71,71) (71,0) (0,0) 

FIG. 3. Comparison of the spin-wave dispersion relations along the symmetry directions of the Brillouin zone, as obtained for 
the itinerant model (full line) and for the Heisenberg model with second (J2/J1 = 0.2) and third (J3/J1 = 0.2) nearest-neighbor 
couplings (broken line). Both dispersions correspond to the same value of Q = 2.556 and are plotted in units of J = it^ /U and 
Ji, respectively. 



where J(q) = Jeff (cosqx + cosqy) and cd{(l)RKKY is given by the expression (|60[). Note, in particular, that if one 
could set uj{q)RKKY = 0, the Heisenberg form ( [5S| ) would result from Eq. (|6l|). 

A final comment on the possible comparison between the spin-wave spectrum we have obtained and|-tiie available 
experimental data is in order. That the classical spin-wave theory (possibly with quantum correctionsE3) can accu- 
rately describe the spin-wave spectrum over the whole Brillouin aonc for the parent compounds of high-temperature 
superconductors (zero doping) is well established at this point .u That when carriers are added (finite doping) in- 
commensurate spin fluctuations occur with rather short coherence length has also been well established in several 
materials.EJ It is interesting to point out, in addition, ^that the observation of well-defined spin waves only close to 
the boundary of the Brillouin zone has been reported,H consistently with what we have obtained by our calculation. 
However, detailed comparison with our form ( |57| ) of the spin-wave dispersion relation with the experimental data may 
not be possible, because our result is valid in the asymptotic limit U » W, where W ^ 8t is the bandwidth, which 
may not be realized for real materials. 



V. CONCLUDING REMARKS 

In this paper, we have studied the spin- wave spectrum for a two-dimensional Hubbard Hamiltonian in the presence 
of an incommensurate spin-spiral phase within the electronic RPA approximation in the broken-symmetry phase. We 
have, in particular, obtained the analytic form of this spectrum at the leading order in the small parameter t/U of the 
Hubbard Hamiltonian. Specifically, it has been possible to obtain the spectrum in a closed form even in the presence 
of an incommensurate structure, owing to the peculiar symmetry which is intrinsic to the spiral phase. In this respect, 
starting from an alternative mean-field configuration different from the spiral one (such as, for instance, a "stripe" 
structure) would have not enabled us to solve for the spin wave spectrum in a closed form. By our approach, we have 
thus pushed the analytic results for the spin wave spectrum of an incommensurate structure as far as possible (apart, 
obviously, from including higher orders in the t/U expansion which could still be done by our approach). 

Alternatively, the spin- wave spectrum could have begn obtained numerically for any value oit/U without performing 
the small t/U expansion, as already reported in Ref.El. In this way, however, it would have been rather difficult (if 
not impossible) to arrive at the functional form ( |57|) for the spin- wave dispersion relation in the strong-coupling limit, 
which has an hybrid form between the dispersion relation for a nearest-neighbor Heisenberg model and that obtained 
within the (long-range) RKKY interaction in the presence of a finite doping. Owing to the itinerant character of the 
system, it is, in fact, the presence of a band of metallic character (crossed by the Fermi level) which generates the 
RKKY magnetic interaction between the localized spins associated with the (filled) lower band. This novel feature 
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constitutes the main result of the present paper. We have concluded accordingly that, even for small doping, the 
itinerant model we have considered results in a dispersion relation w(q) that cannot be effectively represented by 
the Heisenberg model, for which the couplings extend to a few neighbors only. This occurs because the free carriers, 
associated with the itinerant character of our starting Hamiltonian, necessarily introduce long-range RKKY-type 
magnetic interactions among the localized spins. 

A serious concern, which is related to the spin-wave spectrum we have obtained, regards the instability occurring 
about the center of the Brillouin zone, where the spin-wave spectrum becomes purely imaginary. To overcome this 
point, one should possibly start from a more complicated incommensurate mean- field solution other than the spiral 
configuration, with a lower ground-state energy. In this way, however, one would unavoidably not obtain a closed- 
form equation for the spin-wave spectrum, since the incommensurability could in general prevent it. Nonetheless, we 
expect on physical grounds the spin-wave spectrum we have obtained with the spiral pattern to survive inclusion of 
more realistic spin structures, if one considers only the region close to the boundary of the Brillouin zone, for which 
knowledge of the detailed form of the underlying (long-range) spin pattern appears to be less crucial. 
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APPENDIX A: SOLUTION OF THE RPA EQUATIONS FOR THE DYNAMICAL SUSCEPTIBILITIES 
WITH AN INCOMMENSURATE SPIN-SPIRAL GROUND STATE 

In this Appendix, we provide the details of the derivation of the explicit form of the generalized correlation function 
within RPA and of the associated spin- wave dispersion relation, which were reported only schematically in Section 2. 

We begin by showing how Eq. (|l(]| ) of the text can be reduced to a closed-form equation for the correlation function 
X itself. To this end, we note from Eq. that the following relations hold for = and n — z, in the order: 

Xo^u{x, x') = ^ f^a'.fS' L{x+, X f3'; x^-f, x^a) 

a' ,13' 

- ^ E L{x-,x'P'- x+~,x'+a') (Al) 



and 



X^,^{x,x') = E <'..i3'L{x+,x'(3']x++,x'+a') 

a'.p' 



a' ,13' 

By adding and subtracting both sides of the above equations, we obtain 



- ^ x'p'- X+-, x'+a') . (A2) 



Ea' f3' <' 13' a;'/3'; x++, x'+a') = 2 [A'o,,(a;, x') + x')] 

-i Ea'!;3' <'',i3' L{x- , x' /3' ; x+-,x'+a') = 2 [Xo,,{x, x') - x')] 



For jjL — X and n — y we obtain instead: 



J2a',f3' <',f3' x'P'] X++, x'+a') = 2 [X^^Ax, x') + iXy^x, x')] 

Y.a',13' ^a',13' L{x+, x' (3' \ x+ - , x'+ a') = 2 [Xx,y{x, x') - iXy^x, x')] 



(A3) 



(A4) 



with similar relations holding for the non-interacting counterparts of L and X . Entering these relations into Eq. (pi]), 
the closed-form expression (|ll|) results eventually after suitable manipulation. 

Next, we s^cify the form of the non-interacting counterpart A"'"' of X, as defined by the first term on the right-hand 
side of Eq. (|lC|), namely: 
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a, 13 a'. ,13' 



xa) 



lEE</3<',/3' (t U,{x)^l.{x')\){T[i,p,{x'mx)]) 



(A5) 



a,/3 q',/3' 

To this end, we follow Ref.0 and introduce the set of destruction operators along the local spin-quantization 
axes, which identify the spiral pattern of the magnetic ground state within the mean-field approximation. The field 
operator (0) acquires then the form: 



V'a(r) = E E ''^('^ ^ -^^^ 7^(rij)a{ 



(A6) 



BZ 



A'k'R.i 



i i k 

where TZ{U,i) is the spin rotation operator associated with the angles and where the Bloch transformation has been 
introduced as in Section 3, which brings the mean-field Hubbard Hamiltonian into block form for each wave vector k 
belonging to the Brillouin zone {BZ). 

Let Wj_r(k) (with ^ = +, — ) be the matrix which diagonalizes the mean-field Hubbard Hamiltonian at given k, 
such that '[cf. Eq. (|2F 



BZ BZ 

E E 4« C') '^ke = E E ^Ir 7k.. 



k r 



with 

7k,. = 5]<^(k)dk« 

and er(k) given by Eq. (p9|). Upon averaging over the the broken-symmetry ground state, one then obtains: 

(r[dkat)4'e (*')]) -^t^«,.(k)w^;,^^,(k') {T[^^.r{t)llAt') 

where 

= 5k,k' Sr,r' {Qit - t') [1 - fpierim - - t)fF{Am 

Q{t) being the unit step function. Introducing further the tensor T{^i) via the relation 

^7^t(^!,)^a<^m)/3,e =E^M-^("^)'^«' ' 



(A7) 



(A8) 



(A9) 



(AlO) 



(All) 



such that 



\0 -sin6',; 

where 0; = Q • Ri within the spiral-spin pattern we are considering, 111 and approximating 



/I 

cos^i sin^ 

10 

V — sin 9i cos L 



\r~B.,) 0(r-R,) ^ |0(r-R,)|2 5,^. 



(A12) 



(A13) 



for localized atomic (Wannier) orbitals, we obtain the following expression for the frequency Fourier transform of the 
non-interacting correlation function: 
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BZ 
i.'j k,k' 

X Jv^r'(k, k', Oj) . 



Here we have introduced the notation: 



i^/:,,(k,k') ^ E<?(kx«' w^«',.'(k') 



(A14) 
(A15) 



as well as 



J>.r'(k, k', a;) = 



[l-/f(e,(k))]/j.(e,>(kO) 
a; — er(k) + er'(k') + 

[l-/f(6,>(k'))]/F(e.(k)) 
a; — er(k) + er'(k') ~ 



(A16) 



where /i^(e) is the Fermi function. Equation ( [13D of the text is thus recovered. 

We have mentioned in Section 2 that the solution to the integral equation ( pi] ) could be considerably simplified, 
provided the matrix T{U,j) given by Eq. (A12) were preliminary brought to diagonal form by a suitable unitary 
transformation. This transformation reads: 



with a = (0, 1, 2, 3) and where 



B 



E 

/i=(0,a;,y,2) 



/I 

10 
Vo 71 71/ 



(A17) 



(A18) 



In this way, Eq.(A15) is replaced by: 

F,V(k,k') 



^ i?,^F;;,,(k,k') , 

/j=(0,x,y,z) 



and the correlation functions transform according to the rule: 

Xab = E ^'^t^ '^A"' {^^)^b 



The second term on the right-hand side of the integral equation (^T|) then transforms as follows 

BX'-°hXB'^ = B X'-"'^ B^ {B'^y^ e B~^ B XB'^ = X''"hx 
in matrix notation, where 



{B^y' eB-' = 



10 
-1 
0-10 
0-100 



(A19) 



(A20) 



(A21) 



(A22) 



at the same time 



f{n,) = BT{n,)B-^ 



/ 1 \ 

e-*Q-^' 

10 

\ e''^-^' J 



(A23) 
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becomes diagonal, as anticipated. Equation ( p^ of the text then follows. 

There remains to solve Eq. (fiSh of the text explicitly. To this end, we introduce the compact notation 



la 






{a 


= 0) 


-1 


{a 


II 





{a 


= 2) 


+1 


{a 


-3) 



(A24) 



in such a way that (the lattice Fourier transform of) Eq. ( A23| ) reads: 

A(k + 7„Q) (A25) 

A(k) being the lattice Kronecker delta function. In the transformed basis (cf. Eq. ( A20| )), the non-interacting part 
( [l7| ) of the correlation function then takes the form: 



(0)/ 



(A26) 



where we have set 



BZ 



X^^,\ci; iolQ) = ^ E E ^r'A^ - q + 7aQ, k) F,^,,,(k, k - q + 7aQ) 

k r.r' 

X .F„.(k,k-q + 7,Q) . (A27) 
Owing to the wave vector conserving Kronecker delta function in Eq. ( A26[ ) , the integral equation ( p^ becomes: 

Xabiq. q'; = A(q - q' - (7a + 7b)Q) x^°^ (q; '^IQ) 



2U 



4°)(q;L^|Q)A'o;,(q-7aQ,q';^) 

- Xi"^ (q; ^IQ) Xsbiq - ha - 1)Q, q'; c^) 

- xl°^(q;w|Q)i'2f,(q-7aQ,q';^) 

- Xi°^ (q; a;|Q) Xuiq ~ ha + 1)Q, q'; c^) 



(A28) 



Note that the wave vector arguments on the right-hand side of the above expression depend on the index a of the 
matrix element. To avoid this feature, we let q — > q -|- 7aQ everywhere in the above expression and obtain: 

Xabiq + 7aQ, q'; to) = A(q - q' - 7bQ) X^°,\q + 7,Q; c^|Q) 



2U 



xll\q + ^aQ;u;\Q)Xob{<l,q';u;) 

x^ai (q + 7aQ; ^IQ) x^M + Q, q'; ^) 

^l?(q + 7aQ;^|Q)A'26(q,q';c.) 

(q + 7aQ; ^IQ) X^iq. - Q, Q'; c.) 



For any given value of the index h, this relation can then be cast in the equivalent form: 

( x^^crM^^) \ 

[l + X{q,Lu)] 



( '^o6(q,q';'^) ^ 
'^'i6(q- Q,q';w) 

i'2b(q,q';cj) 
Vi'3b(q + Q,q';u;)/ 



xj;;)(q-Q;c.|Q) 



^^"^(q;^IQ) 
V^^"'(q + Q;^IQ)y 



A(q- 



76Q) 



(A29) 



where we have set 

X(q,w) = 2U 
( 



^,$?(q;^IQ) 



^^2^(q;c^lQ) 



^,$;^(q;^IQ) \ 



-Xl°Hq-Q;c^|Q) x(°)(q-Q;c^|Q) X^"' (q - Q; c.|Q) ^^^(q - Q; c^|Q) 



"10 

-X(°)(q;^|Q) 



x(°)(Q;c.|Q) 



xf-^id-MQ) 



V-X(°\q + Q;c^|Q) x(°)(q + Q;c^|Q) x'^^q + Q; u;|Q) x'?^ (Q + Q; c^|Q) / 



(A30) 
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Solving for Xab in Eq. ( [A29| ), we obtain eventually: 

Xabiq + 7aQ, q'; ^) - ^ [1 + X(q, u;)]-^, 

a' 

X x(?)(q + 7,,Q;L^|Q)A(q-q'-7bQ) 

which coincides with Eq. ( [l9| ) of the text. 



(A31) 



APPENDIX B: DETAILS OF THE t/U EXPANSION 

In this Appendix, we provide details of the t/U expansion of the mean- field parameters and of the matrix elements 
of the correlation function which are necessary to obtain the spin-wave dispersion relation. 

We begin by considering the leading (t — 0) term of the expansion (^5|) for the mean-field (band) eigenvalues, 
corresponding to completely flat bands, which is given by [cf. Eq. ([2S|)] 

(Bl) 



e,(k) = U (e(") - /.(°)) - {/(mi + {-lym^) - A^o 



with r = 1, 2 and where /ip = [/ /i'"-* is the chemical potential at the lowest order in t/U. Since at this order the band 
eigenvalues do not depend on the wave vector, there is no value of the chemical potential fj,^^"' satisfying equation ( |3^ ) 
at zero temperature and noninteger doping S. More precisely, the doping jumps from 5 = —1 when ^^^'^ < e^i\ to 
6 = when ef ^ < /z''^-' < and to S = +1 when /^^°) > In order to consider a continuously doped system, it is 
therefore necessary to include at the outset the next-to-leading order in the expansion ( ^5|) of the band eigenvalues, 
which introduces a band dispersion. 

For definiteness, we shall consider the case 5 > from now on. (The case S < can be recovered by exploiting 
particle-hole symmetry.) In the case d > 0, the lowest band is always filled, and we will consistently use (the zero- 
temperature value) /^(ei(k)) = 1 throughout. On the other hand, the zeroth order chemical potential for < i5 < 1 

must be chosen as /i^*^-* = e2°\ and the Fermi function in the upper band can be expanded as 



/F(e2(k)) 
<5(.W(k). 



(1) 



e(4^'(k) 

:) (.r(k)-M(^))+(| 



(B2) 



(.f(k)-M(3)) 



--5'(.«(k)-^W 



i) (er^(k)-M^^))^ + 



where we have considered the zero-temperature limit of the Fermi function and of its derivatives, and introduced the 
step function 8 as well as the Dirac S{x) function (not to be confused with doping). In this way, Eq. (^) can be 
expressed in powers of t/U, yielding 



BZ 



(1) 



1 



^Ee(4"M 

k 

BZ 

-1^5(e«(k)-^W)(4^)(k)-M(^)) = 0, 



k 

BZ 



^ Y: [s' (e«(k) - M^^') (4^'(k) - + 2 S (4^^(k) - ^(^)) (4^^(k) - M^^^) 



= 0. 



Similarly, the small t/U expansion (36) of Eq. ( |33| ) yields [cf. also Eqs. (37)-(p9| 




(B3) 



(B4) 



(l-/^^(e2(k)))+0 - 
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where the Fermi function should be again expanded in the form (B2). One obtains eventually: 



,(0) 



.(1) 



.(2) 



BZ 

6 (4"(k) -...,), 



k 

BZ 



k 

BZ 



A/" 51 



ro(k) 



BZ 



(0) 



e(.W(k)-;.«" 
1^£<5' (e«(k)-^W)(4^)(k)-M(2))^ 

k 

BZ 

^5:<5(.«(k)-^w) (.f(k)-M(3)) 



k 

BZ 



ro(k) 



2m. 



(0) 



e(.W(k)-M«) , 



(B5) 



where in the last step we have used the last of Eqs. ([B^). 

There remains to expand in powers of t/U the last the self-consistency equation (^^. For the diagonal spin-spiral 
solution we obtain: 



^ BZ 

= sin(g/2)-^^ {coskx + cos ky) (1 -I- /F(e2(k))) 



cos(Q/2)- 



(B6) 



1 1 



BZ 



J2 Toik) (sinfc, + sinfcy) (1 - /F(e2(k))) 



where the expansion (B2) has still to be inserted. 

From Eq. (p9[), the argument of the 8 function is T'e(k) +7712^-* — /i*^^^ (with Te(k) = 2 cos((5/2) (cos/c^; +cosA;j,) for 
the diagonal spiral solution), which can take positive as w ell a s neg ativ e values. 

At the various orders in t/U we then obtain from Eqs. ( p33|) and (|BE|): 
(i) order {t/U f 



BZ 



1 - e (Te(k) +m 



(1) - „(i) 



.(0) 



BZ 



(B7) 



(ii) order {t/Uy 



BZ 

1 5: S (T,(k) + - M^^O (4''('^) - ^^'^^ = « ' 



.(1) 



(B8) 



(iii) order {t/Uf 



J2) 



^ BZ 



To(k) 



(0) 



e(re(k) + r4'^-^(i) 



(B9) 
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and we don't need to evaluate {e'^'' (k) — fi'^^^) at the order we are considering. It is convenient to solve these equations 



,(3) 



order-by-order for the two variables /i*-*-* 



and 



(z = 0,l,2). 



Using at this point the method developed in Appendix C for performing the k summation over the relevant portions 
of the Brillouin zone, we obtain eventually the following results for the mean-field parameters at the leading orders 
in S: 



4°^ - 
,(1) - 



m 



(2) 



and 





(0) 




(1) 

711,2 


^ ^(^) - 


(2) 
7172 



(l-<5)/2 


-4 sin2(Q/2)+0(5), 



mi = (l + (5)/2 
4cos(Q/2)(l-7r(5 + 0((52)) 
Q + 0{S) . 



(BIO) 



(Bll) 



There remains to determine the magnitude Q of the characteristic wave vector from the self-consistency condition 



(B6), where for the Fermi function we use the expansion (B2). Equation (B6) then becomes: 

BZ 



sin(Q/2) (CO" 

H- COS Ky 

k 



(1) ( 

7172 — 



2-e re(k) 



,(2) 



,(1) 



(1) 



1 1 



,(1) 



-) cos(Q/2) 



/I) 



2m^ 
= . 

In particular, at the lowest order in t/U we obtain: 

sin(Q/2) {2S + 0{S^)) 







(B12) 



(B13) 



For fixed S and large U/t, one thus has only the ferromagnetic solution Q ~ 0. However, if one allows S to be of the 
order t/U, Eq.(B13) has to be considered together with the next term in t/U, and we obtain instead: 







(B14) 



sin(Q/2) {2 6 + 0(6'^)) + 4 sin(Q/2) cos(Q/2) (1 + 0{6)) 
which yields the two solutions 

sin((5/2) — (ferromagnet) 
cos((5/2) = —SU/2t (diagonal spiral) , 

with the spiral solution being energetically favored over the ferromagnetic solution whenever it exists. From the above 
equation, it is clear that the spiral solution exists for any 6 < 2t/U, which is consistent with our assumption that the 
doping parameter 6 is at most of order t/U in the spiral phase, thus justifying the expansion of Appendix C. The 
transition to the ferromagnetic state is second order, as the incommensurability vector Q decreases continuously from 
Q — TT to Q = 2n with increasing S. 

We pass finally to consider the t/U expansion of the matrix elements of the correlation function. To begin with, we 
use the expressions (|37|)-(|39|) for the eigenvectors of the mean- field Hamiltonian and obtain from the definitions (|22" 
and (H): 



BZ 



"(q'^) = ^E 



l-(-) (T2(k)+T„2(k-q)) 



J^2,i(k, k-q,w) 



[r2(k).Fi4(k, k - q, c^) + T^{k - q).F2,2(k, k - q, u;)] 



o({t/uy 



(B15) 
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+ [J^i,2(k, k- q,cj) + J^2,i(k,k- q,tj) 

- ^i,i(k,k-q,c^)-^2,2(k,k-q,c^)] ]+0({t/Uf 



(B16) 

for the matrix elements entering the final expression (|5^) of the spin-wave dispersion relation. For the remaining 
matrix elements we obtain instead: 



X^-'iq, = ^ 2 I (^i^ (ro(k) - T„(k - q))^ 
X [J^i.2(k,k-q,w)+J^2,i(k,k-q,cj)] 

2 



l-(^) (T,(k)-T„(k-q)f 



X [J^i.i(k,k-q,tj)+J=-2,2(k,k-q,u;)] 



X^'-{q,u;) = ^V2 |(^^) ^ (T^Ck - q) - T<,(k)) ^2,i(k, k - q, 

[r„(k) Ti.iik, k - q, - To{k - q) ^2,2(k, k - q, iu)] 



o 



k 

t \ 1 



U J 1-5 

o[{t/uf) , 



A'^(q,c.) = 



X [J^2,i(k,k-q,cj)- J^i^2(k,k-q,w) 

2 



+ 



l-(-) (T„^(k)+r„2(k-q)) 



[^i,i(k, k - q, c^) - .F2^2(k, k - q, w)] \+0{ (t/U) 



(B17) 



(B18) 



(B19) 



A'o+'^(q,c.) = 2_Y^,^!^(^l^l_L_ (r,(k-q)+T,(k)) .F2,i(k,k-q,c.) 

[~r„(k) ^1,1 (k, k - q, c^) - To(k - q) ^2,2(k, k - q, to)] 



k 

t \ 1 



U J 1-5 

o{it/uf) , 



(B20) 



r (q. -) = ^ E { (^) ' (k) + 7^o(k - q))^ 

X [J^i,2(k,k-q,tj)+J^2.l(k,k-q,t^)] 
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l-(^) (T„(k) + r„(k-q)f 



X [J^i,i(k,k-q,cj)+JC-2.2(k,kk-q,tj)] 



(B21) 



oiit/uy 



Note that in the above expressions J-"i_i(k, k — q, cj) =0 when the doping parameter S > 0. In addition, JFi^2 can be 
obtained from ^2,1 through the symmetry condition (44). Making use of the expansion (^5|) for the eigenvalues of the 
mean-field Hamiltonian, we write further 



as well as 



J-2.2(k,k - q) = 1 [/;^(e2(k - q)) - (e2(k))] — \ 

t re(k - q) - Te(k) 

+ 0{1/U) 



k - q) = [/f (£2(k)) - 1] ^ { (1 „ ] „ ^) 



[T,(k)-Te(k-q)] 



Uj (l-^-J))2 

T2(k)-r„2(k-q)-2m(^) 



1 



(l-5-J))3 



(re(k)-re(k-q))' 



(B22) 



(B23) 



with uj = Lo/U ^ where the expansion (]B2|) of the Fermi function has still to be inserted. 



APPENDIX C: SUMS OVER THE TWO-DIMENSIONAL BZ IN THE SMALL 5 LIMIT 



In this Appendix, we develop a method suitable for performing the k summation over the relevant portions of the 
Brillouin zone in powers of the doping parameter 5. This is, in turn, justified by the fact that 5 in the spiral phase is 
of order t/U , as shown in Appendix B. 

The typical integral to be evaluated is of the form: 



BZ 



(CI) 



where Telk) — 2 cos((5/2)(cos fc^, + cosfc^) for the diagonal spiral solution we are considering and 7 is a parameter 
(which depends on the chemical potential //'■^•' at the order O \{t/U)^\) that in the following we shall simply call /i 
for simplicity. Introducing the notation such that 7 = 4cos((5/2)(l — /^/4) for any given Q value, recalling the 
definition of the 6 function, and considering that cos(Q/2) < in our solution, we set further: 



E 



ff(k). 



(C2) 



cos k^+cos fcy<2— /^(/2 



Quite generally, we can introduce a polar representation for the two-dimensional wave vector k and determine for 
each value of the polar angle </> the magnitude fc(0), such that the equality 



f I 

cos (fc((/)) cos (/>) -f cos sin 0) = 2 — 



(C3) 



is satisfied. In this way we rewrite the integral (C2) as follows: 



If 1 

'^^'^'^ ^ 4^ y^^*2;rffcj,g(fc:r,fcy) - — 2 



fe(0) 



dk k g(k, (p) , 



(C4) 
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with a slight (albeit harmless) abuse of notation for the function g. 

In practice, when is small compared to unity (as it is the case for small 5) it is convenient to determine k{(j)) in 
powers of /^t by expanding cos (fc(</)) cos (p) + cos (A;(0) sin </>) in a Taylor series about fc = 0. For instance, at the lowest 
order we obtain from Eq. (|C3| ): 

fc(0)2 + O {k{c^f) = (C5) 
which gives fc(0)^ = + O (/^)- At the next significant order we obtain instead: 

(sin'' (j) + cof** -^^ 



12 



■fc(0)4 + 0{k{c^f) = 



which gives 



k{^Y 



.fu + 



(sin^ 



cos 



12 



(C6) 



(C7) 



Consider, for instance, the first of Eqs. (B7), which is of the form (CI) with gCk) = 1 and 7 = 7712^'' — /Lt^^-'. Equation 
(C4 ) now becomes: 



Introducing the notation as above, for small values of 6(> 0) we obtain: 



where use has been made of Eq. ( |C7| ). Inverting this relation we obtain eventually: 

= AttS - TT^S^ + O (,53) , 



(C8) 



(C9) 



(CIO) 



which confirms the fact that = 0{S). 

In general, the case of a k-d ependent g can be treated by expanding g(k) about |k| = 0, since the last integral on 
the rignt-hand side of Eq. ( p4|) is restricted to |k| < k{cf)) — 0{S^^'^), and by retaining the relevant orders in i5 in the 
final expression. The integral over the whole BZ (i.e., the first term on the rignt-hand side of Eq. (p4D), on the other 
hand, can usually be performed analytically. 

The Brillouin zone sums involving the Dirac delta functions or its derivatives can be further evaluated according to 
the following device. Writing 



5(Te(k) + 7) =|-e(r,(k) + 7) 



d 



cos(0/2) a/p 



e(Te(k) + 7) 



we obtain for the Brillouin zone sum 

BZ 



;^^.g(k)5(") (r,(k) 



,(1) 



(1) 



(Cll) 



(C12) 



(cos(Q/2))("+i) a(///) )("+!) 
1 ^("+1 



T-rY 9(k) e (Te(k) 



,(1) 



(cos(Q/2))(»+i) a(/(i))(n+i) 47r2 ''Jo 



27r 



fe(0) 



dkk g{k, (j)) 



since the first term on the right-hand side of Eq. (C4) does not depend on /p. 
As an example, we evaluate the expression 



23 



BZ 



which appears in the first of Eqs. (B8). According to Eq. (Cll), we first evaluate the expression (cf. also Eq. (C4)): 



BZ 



(C13) 



4 sin2(Q/2) - -!jsin2(g/2) / d0 / dfc fc (sin(fccos0) + sin(fcsin0))' 



^ JQ JO 

Expanding the factor within parentheses in the last integral in powers of fc, as explained below Eq. ( |C1C| ), we obtain: 

BZ 



;^^r„2(k)e(r,(k) + mW-^(i: 



1 /-Stt pk{(t)) 

= 4 sin^ (Q/2) - — sin^ {Q/2) dcj) dkk {k^ (cos </> + sin + O (/c^) ) 



"'0 



3in2(Q/2) {A-8ttS^ + 0(6^)) 



Using Eq. (Cll) we obtain eventually: 

BZ 

^^r„2(k)<5(r,(k) + ™«-;.«) 



(C14) 



(C15) 



1 d 1 

cos(Q/2) a/^ 77 ^ 



BZ 

^r2(k)e(r,(k) + mW-A*(i: 



__sin^l ,^ sin^(Q/2) 

" cos(g/2) TT + ^ ^^-^^'^ " cos(0/2) + ^^'^ ^ 

With these prescriptions, the results ( |B10| ), ( Bll ), and ( |B12 ) are readily obtained. 



It was shown in Section 3 that, as far as the matrix elements of the correlation function are concerned, only the 
explicit expressions of d(q, w) and e(q, w) are required at the order in t/U we are considering in this paper. Writing 



and 



= /(q, Lo) - J(q, Lo) - K{q, uj) + M (q, lo) 
e{q,uj) = L{q,uj) + L{-q,-uj) + N{q,u}) 



(C16) 



(C17) 



where the quantities I, J, ■ ■ ■ are specified below, we obtain for the sums over the wave vector k using the method 
described in this Appendix: 



U 



BZ 



= JfYl •^2i(k,k - q,( 



-l + S 



+ ( — ) 2S cos(Q/2) (cosg^, + cos qy) 



(C18) 



— 4 cos- 
U 



'■{Q/2){cosq.,+cosqy) + o(^it/uf^ , 



2 BZ 



J(q,c.)^^(i) J2Toi^)^2iiK^-ci,u;) 



-4 I — I sm 

U 



m\Q/2) + o[it/uf) , 



(C19) 



24 



if(q,c^) 



U t 



2 BZ 



T2(k-q)^2l(k,k-q,L^) 



= -4 



^ysm\Q/2)+o[it/uf) , 



(C20) 



L{q,u) 



U ( t 



Af \U 



2 BZ 



T„(k)ro(k-q)^2l(k,k-q, 



-2 ( - ) sin2(Q/2) (cosg:, + cosgj,) + O ( (t/C/)' 



(C21) 



2 BZ 



M(q,c.)^-^ (Ij ^r2(k-q)^22(k,k-q,c.) 

k 

-26 il) sin^jQ/2)(sing. + sing,)^^^^^/^^^^3 



[/ y cos(Q/2) (cos(7j. + cos^j, — 2) 



(C22) 



N{q,u;)^^ (^) ^To(k)To(k-q)^22(k,k-q,o.) 



A/" \U 



2 BZ 



where these results hold for |q| :s> k p as discus sed i n Section 3-B. 

For the relevant matrix elements ( C16| ) and ( C17 ) we obtain eventually : 



d(q, u) ~ —I — uj - 



^\ A cos^ {Q / 2) {cos q^,+ cos qy- 2) +^(^\ sin2(g/2) 



2 I —j 5 cos(Q/2) (cos 9a; + cosgy - 2) 

U ) cos(Q/2) (cos^x + cos9j^ - 2) ' 



and 



(C23) 



(C24) 



e(q,w) 



-4 ( - ) sin2(Q/2) (cosQa; +cosgy) + 0[{t/Uy 



(C25) 



APPENDIX D: HEISENBERG SPIN WAVES FOR A SPIRAL CONFIGURATION 

In this Appendix, we give a simplified derivation of the spin-wave dispersion relation for a two-dimensional Heisen- 
berg antiferromagnet iiJ-tJi^ presence of a spiral incommensurate magnetic ground state. Although the form of this 
spectrum is well known," 'irH it is worthed to give here a compact derivation in terms of the set of local spin quantization 
axes utilized in Section 3 and in Appendix A. 

We consider the Heisenberg Hamiltonian 

^ = ^E-^-^S.-S, (Dl) 

for spin 1/2, where the sum over the lattice sites i and j extends in principle over distant neighbors. We transform 
next the spin locally as follows: 
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where T^"^^ represents the 3x3 matrix given by the spin part of (|A12|). The Haniiltonian (Dl) then becomes: 



(D2) 



(D3) 



where T^i^^ contains the difference Qj — 9i in the place of 9i of Eq. ( A12 ). 
Writing further 

= \ {St + S-) , Sf = ^ [St Si) , (D4) 
and carrying out the standard Holstein-Primakoff transformation 

Sf ^ V2Saf 

5f — alui — S = Hi — S (D5) 
at the leading order away from perfect alignment (and for large values of S), the Hamiltonian ( |D3D becomes: 

52 



iJ = — ^ Jij cos{dj - 6'i) - S* ^ Ji J cos(6'j - 9i) 

+ f E -^'j + «]«*)[cos(f?, - e,) + 1] 

+ (a|a] + aiaj)[cos{dj - 6,) - 1]} ■ 



We introduce at this point the lattice Fourier transform, specify the choice 0i = Q • Ri, and write: 

BZ 

2 



^ BZ 

^ a^hj cos Q • (Rj - R,;) = x E a(k)6(-k) ( J(k - Q) + J(k + Q)) , 



where and bj is an arbitrary pair of operators and where 

J(k) = exp(-ik • R) J(R) 



R 



With the notation 



and 



o S*/ , J(k + Q) + J(k- 

a(k) = -5 J(Q) + - f J(k) + ^ ^' ^ ^ 



/3(k).^f(j(k)- ^(^+Q) + ^(^-Q) 



(D6) 



(D7) 



(D8) 



(D9) 



(DIO) 



for any given value of Q, the Hamiltonian (D6) reduces to the form: 

BZ 



H = 



AfS^ 



J(Q) + ["(k)a^(k)a(k) 



+ ^(at(k)at(-k) +a(k)a(-k)) 



(Dll) 



The standard Bogoliubov transformation can be used at this point to diagonalize the Hamiltonian ( Dll ). Writing 

a(k) = u(k)6(k) + w(k)6^(-k) (D12) 
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where M(k) = cosh'0(k) and w(k) = sinh'!/'(k) with 



tanh2?/.(k) = -ffl = 7(k) , 



a(k) 



the Hamiltonian ( Dll ) becomes eventuaUy 



BZ 



H ^ Eo + El +Y, e{k)b^k)b{k) 



(D13) 



(D14) 



In this expression: 



(D15) 



^1 = - E («(k)^^'(k) + u(k)z;(k)/3(k)) , 



(D16) 



and 



£(k) = v/a'(k) - /32(k) 



5 W(J(Q)- J(k)) J(Q) 



J(k + Q) + J(k-Q) 



(D17) 



The spectrum ( pl7 ) coincides (apa rt fo r a different normahzation of J) with that given in Ref.0 (cf. in particular 
Eq. (67) therein). Note from Eq. ( D17 ) the presence of a Goldstone mode for k = and k = ±Q (owing to the 
symmetry J(— k) = J(k)). 

In the classical limit (that is, for large values of S), minimization of the ground-state energy is equivalent to finding 
the minimum of J(Q). In this way, the appropriate value of Q is determined. For instance, when second (J2) and 
third (J3) nearest- neighbor couplings are included besides the nearest- neighbor antiferromagnetic coupling (Ji), J(Q) 
reads: 



J(Q) = Ji (cosQa; -l-cosQj^) -I- 2J2 cos Qa; cos + J3 (cos2Qa; + cos 2(5 j^) 



(D18) 



(with the wave vectors measured in units of the inverse of the lattice spacing). In this case, the diagonal spiral 
configuration is stable for J3 > Ji/4 — J2/2 and J3 > J2/2, implying that J3 must be nonvanishing for the spiral 
phase to be stable.E3 The corresponding spin- wave spectrum ( |D17| ) for J3 = J2 = Ji/5 and S = 1/2 is shown in Fig. 4 
(full line), where the antiferromagnetic spectrum with Q — (tt, tt) and J3 = J2 = (broken line) is also shown for 
comparison. 



APPENDIX E: RKKY SPIN WAVES FOR A SPIRAL CONFIGURATION 



In this Appendix, we show how the part of the dispersion relation (|5^) of the text containing the trigonometric 
functions in the denominator originates from the Rudermap-Kittel-Kasuya-Yosida (RKKY) interaction between two 
localized spins and mediated by the conduction electrons .r^W ^1 The long-range nature of this interaction is readily 



evidenced by expanding formally the denominator in Eq. (57) as a a power series in cosg^, and cos qy, so that an 
infinite number of Heisenberg-like terms appaars. 

We follow here the treatment by Mattisj23 and consider the interaction Hamiltonian between the localized spin 
operators (associated with the valence (filled) band) and the itinerant spin operators Sc(Ri) (associated with the 
conduction band), where the sufRx i specifies the lattice site. We write accordingly: 

i 

where / is an exchange integral. 

For a generic spin-1/2 operator at site i associated with the conduction electrons wc write: 
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2.0 

co(q) 
Jl 

1.0 




(0,0) (7t,7t) (7t,0) (0,0) 

FIG. 4. Spin-wave spectrum for the Heisenberg Hamiltonian with a diagonal-spiral spin configuration and couplings up to 
third nearest neighbors specified in the text (full line) . The antiferromagnetic spectrum with nearest-neighbor coupling only is 
also shown for comparison (broken line). 



where fj, = {+,—,z) and di^ are the destruction operators along the local spin-qua ntiz ation axes. Combining these 
opera tors with the eigenvectors of the mean-field Hubbard Hamiltonian [cf. Eq. (|A8|)] and recalling the notation 
( A15| ), we cast Eq. ( |E^ ) in the form: 



- ^ E E ^k+q,. Pry (k + q, k) 7k,.' . (E3) 



k,q r,r' 



The restriction to the conduction band implies that r = r' = 2 in Eq. (E3) when 5 > 0. In this case the Hamiltonian 
(El) becomes: 



^ EE («r)' (e ^-^''•"'7U,2^2^2(k+q,k)7k,. J 

i H \k,q J 



iJcxc - > . > . (Sf)' I > . e--!-' 7^0.2 ^2"2(k + q, k) 7k,2 | • (E4) 

Applying at this point the standard procedure^ to obtain the energy shift of second-order in the Hamiltonian (^), 
the following RKKY-type effective Hamiltonian results: 

i?™ = EE(^n'4''^j (E5) 

ij iiu 

with the notation 

J^^- ^ \- \- ,-^q.(R.-R,) j^2"2(k + q,k)Fj-;(k-fq,k) 

-^^1 - 4^2 2^ 2^ e2(k) - e2(k + q) ' ^ ' 

k<kF |k+q|>fcF ' '^^ ^' 



Note that the expression (E6) does not vanish even when sites i and j are far apart. Note also that the energy 
denominator on the right-hand side never vanish by construction. 

There remains to obtain the spin- wave dispersion associated with the spin Hamiltonian ( |E5| ). To this end, we 
introduce in Eq. (^) the usual Holstein-Primakoff transform ation for the spin operators and truncate the expansion 
in the bosonic operators a and to quadratic order [cf. Eq. ( |D5| )], to obtain: 

i/«^^^ = £;„ + 5 E [(^++(q) + ^"~(-q))«^(q)«(q) 

|q|»feF 

+ J+-(q)at(q)at(-q) + J-+(-q) a(q) a(-q)] (E7) 
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where now 



|q|>fci? |q|>fcF 



(E8) 



with the restriction |q| ^ kp replacing the weaker condition |k + q| > kp (cf. Section 3-B). Note that in Eq. (E7) we 
have introduced the notation 



(E9) 



and used the identity J2i JiT — for any given j. 



The quadratic Hamiltonian (E7) is then diagonahzed by the standard BogoHubov transformation (cf. Appendix 
D). The resulting spin- wave spectrum has the form: 



UJ 



RKKY t 



(q) = S (j++(q) + J-(_q))2 _ 4 J+- (q) J-+ (-q) 



1/2 



(ElO) 



To proceed further, we need to specify the quantities Fj^, -F22' ^2 entering the definition (EE), at the relevant 
order in the small parameter t/U. Using Eqs. (pa)-(pa) we obtain for i5 > 0: 



F2+2(k + q,k)=zV2 (-) To(k + q) + O {{t/Uf'' 



F22(k + q,k) 



iV2 (^1^ To(k) +0((i/t/)2) 



as well as 



e2(k) - e2(k + q) = t {T,[\<i) - T„(k + q)) + O {[t/Uf) , 



(Ell) 



(E12) 



where we have approximated 27712*^' = 1 — (5 = 1 



The method developed in Appendix C can be used at this point to perform the sum over k in Eq. ( |E9| ). The result 



is: 



^ ^2 i j Sin^(Q/2) (sing^ +singy)^ 

U"^ cos((5/2) cos (7a; + cos qy — 2 ' 



(E13) 



J-(q) - -/2_0(J2)^ 



(E14) 



J+-(q) = J-+(q) - I^^O{6'). 



At the lowest order, only J^^ contributes to Eq. (ElO), which then reduces to: 



'"^^'^fa) ^ cj2 t c sin^(Q/2) (sing^ +singy)^ 

C/2 |cos(g/2)| (2- cosq:r - cosg^) 



(E15) 



(E16) 



The self-consistency equation (^) relating Q to 5 and t/U for the diagonal spiral configuration can eventually be 
used, to yield (for S" = 1/2): 



.RKKY , 



q) = JcS 



(singa; -t- singy)^ 
2 — cos qx — cos qy 



(E17) 



with Jeff given by Eq. ( ^8[ ) of the text and where we have set / ~ 2U. Note that Eq. ( E17 ) coincides with Eq. (|60| ) of 
t he te xt, which was obtained by setting "by hand" T21 = T\2 = in the full calculation. Note also that the expression 
(E17) involves transverse spin components only, akin the general expression (53) of the text. 
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